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We perform a quantum mechanical analysis of superconducting resonators subject to dielectric 
loss arising from charged two-level systems. We present numerical and analytical descriptions of 
the dynamics of energy decay from the resonator within the Jaynes-Cummings model. Our analysis 
allows us to distinguish the strong and weak coupling regimes of the model and to describe within 
each regime cases where the two-level system is unsaturated or saturated. We find that the quantum 
theory agrees with the classical model for weak coupling. However, for strong coupling the quantum 
theory predicts lower loss than the classical theory in the unsaturated regime. Also, in contrast to 
the classical theory, the photon number at which saturation occurs in the strong coupling quantum 
theory is independent of the coupling between the resonator and the two-level system. 



I. INTRODUCTION 

Noise is a central issue in the implementation of quan- 
tum computation using superconducting circuitsii^. Re- 
searchers have focused on two-level charge fluctuators^ 
as candidate sources of energy- and coherence-loss in su- 
perconducting qubits^"— . Superconducting resonators, 
which are sensitive probes of bulk two-level system (TLS) 
loss, are also useful for dispersive qubit readout using cir- 
cuit QEDi^, single artificial-atom lasingi^ and single pho- 
ton detection^. Experiments focusing on dielectric loss 
in superconducting resonators have previously employed 
a classical description of the resonato r^^i^^ . However, 
many of these experiments involve low microwave-photon 
numbers and properly require a full quantum mechanical 
analysis for their interpretatio n^^i^^ . 

In this article we quantize the classical theory of loss 
due to two-level systems (TLSs) to obtain the famil- 
iar Jaynes-Cummings model ( JCM) of quantum optics^. 
Although the JCM has been studied thoroughly in both 
atomic and solid-state physics, to the best of our knowl- 
edge it has not been used to explore the loss regimes 
presented her o^°i^^ , nor has it specifically been used to 
evaluate loss tangents. (In a similar way, there has 
been recent renewed interest in the JCM arising from 
nanofabricated solid state systems that are beginning to 
broach parameter regimes inaccessible to atomic physics 



experimenti 
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In this paper we use the JCM to analyse a simplified 
set-up designed to capture the essence of resonator loss : 
a harmonic oscillator that can release quanta to a zero- 
temperature thermal bath only through a resonant TLS. 
We neglect the fact that a resonator may interact with 
many TLSs. Working with this minimal model allows us 
to establish the dynamics of decay and to demarcate var- 
ious parameter regimes. This model can serve as a start- 
ing point for studies of quantum mechanical loss in more 
complex systems. It also allows us to generate transpar- 
ent analytic expressions for the resonator photon number, 
TLS excitation and coherence, as well as system correla- 
tions, which are of interest not only in the context of res- 



onator loss but also for control and manipulation of the 
full system^ii^. Our analysis shows agreement between 
the quantum and classical theories for weak coupling be- 
tween the resonator and the TLS. However in the case 
of strong coupling, we find the quantum theory predicts 
a significantly lower loss than the classical theory when 
the TLS is unsaturated. Also, in contrast to the classical 
case, the photon number at which the loss saturates does 
not depend on the resonator-TLS coupling. 

The remainder of the paper is arranged as follows. In 
Section |IT] we summarize the classical model of loss and 
proceed to quantize it, arriving at the Jaynes-Cummings 
model; in Section IIIII we describe two numerical ap- 
proaches towards solving the problem, using the master 
equation and the Maxwell-Bloch equations respectively; 
in Section IIVI we present approximate analytical solu- 
tions to the Maxwell-Bloch equations which allow us to 
identify loss regimes exhibiting qualitatively different be- 
havior. Section |V] supplies a conclusion. 



II. CLASSICAL AND QUANTUM MODELS OF 
LOSS 



In this section we briefly restate the TLS model that 
treats the superconducting resonator classically and pro- 
ceed to quantize it to obtain the JCM. This exercise al- 
lows us to directly relate the parameters of the two mod- 
els. 



A. Classical model 

The classical model consists of a defect with charge 
q that can tunnel between two spatially distinct sites 
separated by a distance / inside the dielectric that is 
permeated by the electric held of the superconducting 
resonator—. The onsite energies of the two states \L) 
and \R) are ±A/2 respectively and the height of the tun- 
neling barrier is Ag. The defect is driven by an electric 
held of magnitude F{t) which is directed at an angle 9 
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with respect to a line joining the two charge sites. The 
classical Hamiltonian of this system may be written as 
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E p cos 9F{t) A 
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(1) 

where = A"^ + Aq,p = ql and the Pauli matrices <Tz 
and ax have been defined in the energy basis 



1+) = cos||L)+sin||i?) 
h> = sin||L)-cos||i?), 



(2) 



with tana = Aq/A. The last term in Eq. ([T]) denotes 
the field energy contained in the resonator 



U 



eF\t) 



B\t) 



(3) 



where e and /i are the electric permittivity and mag- 
netic susceptibility respectively and B{t) is the magnetic 
field. The energy U is usually not included in the clas- 
sical model since the fields are not dynamical variables 
in that case; however it is relevant to the quantization of 
the problem below. 



B. Quantum model 

To quantize the electromagnetic field in the Hamilto- 
nian of Eq. (HI) we use the prescriptions^^ 



U TiLoa'^a, F{t) F^ [a{t) + a'^{t)] 
where to is the frequency of the resonator, 



(4) 



which is of the form of the standard JCM, Hjcm 
hg{a(T^ + a^(T_), with a coupling constant given by 



9 



pcosOF^Ao 
2hE ■ 



(8) 



In the remainder of the paper we will persist in using g 
for the sake of avoiding lengthy expressions. Below we 
will consider the evolution of the TLS-resonator system 
when it is connected to reservoirs that cause relaxation 
and dephasing. 



III. EQUATIONS OF MOTION : NUMERICAL 
TREATMENT 

In this section we will consider the full quantum me- 
chanical treatment of the problem numerically via solu- 
tion of the density matrix, and also present an approxi- 
mate but simpler numerical approach using the equations 
for the expectation values of the relevant physical quan- 
tities. 



A. Master equation 

The presence of dissipation and dephasing in our prob- 
lem can be accounted for by using the standard mas- 
ter equation approach for the JCM, which yields, in the 
Born-Markov approximation, an equation of motion for 
the density matrix p of the TLS-resonator systero^i, 

p = -ig[a^(j^ + aa+,p] + -^{(Jzpaz - p) 

+ 7-^(2cr_pcr+ -cr+cr_p-p(7+cr_). (9) 
11 1 
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is the 'electric field per photon', V is the effective mode 
volume of the resonator and 0(0^) is an annihilation (cre- 
ation) operator obeying the bosonic commmutation rule 
[a, a^] = 1. We thus obtain the quantum Hamiltonian, 



He 



-+pcos0F^ — (a + flt) 



where a± are the TLS raising and lowering operators 
respectively. Assuming the resonance condition Tilj = 
E, transforming to the interaction picture with respect 



to Hq = Tiuiia 



^) and making the rotating wave 



approximation, we arrive at 

H'q ^p cos 6 F^^{aa+ 



(5) The first term signifies unitary evolution, the second cor- 
responds to coupling to a reservoir that causes pure de- 
phasing of the TLS at the rate l/T^ and the last term 
denotes a coupling with a zero-temperature reservoir into 
which the TLS can relax at the rate 1 /Ti . 

Our approach will be to begin with the superconduct- 
ing resonator in a coherent state |a) (with average photon 
number jap = (n(0))) and the TLS in its ground state, 
and study the rate at which quanta are lost from the 
resonator. For simplicity, we have not included an exter- 
nal drive for the resonator. We have also assumed that 
the resonator suffers no intrinsic loss since we wish to 
study loss via the TLSsi^i^. The typical frequencies of 
resonator operation correspond to energies much higher 
than available in a cryogenic environment {Ulo ^ ksT), 
justifying our assumption of a zero-temperature reservoir 
above. Other regimes of the damped TLS-oscillator sys- 
tem have been addressed in, for example Ref. 

The numerical solutions to Eq. ([9]) wiU be discussed in 
the context of the analytic results presented in Section 
(J) EV] below. 



(6) 
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B. Maxwell-Bloch equations 

While Eq. ^ allows us to obtain a fully quantum 
mechanical solution to the problem, it can involve the 
population and coherence dynamics of a large number of 
states, especially for high initial photon numbers in the 
resonator. In this case it is useful to have a less intensive, 
and only slightly less rigorous, numerical approach to the 
problem. This route is provided by the Maxwell-Bloch 
equationsS^ which follow from Eq. ((9)) , 



ig{a) 



To 



(10) 
(11) 



= -ig ((aa+) - {a^a^)) - — (12) 



where <J++ = (1 + (Tz)/2 is the population operator for 
the upper TLS level |+), and 
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We can also find from Eq. © the equation for the dy- 
namics of the average photon number. 



{n{t)) = -(ata) 



(ata_)) 



Combining Eqs. ([T2l) and (fT4|) we find, 

^[(n(t)) + (aH.+)] = -^(a++), 



(14) 



(15) 



which represents the conservation of energy in the sys- 
tem. The LHS of Eq. ((T5|) denotes the rate of change of 
the the combined TLS-oscillator energy, while the RHS 
signifies the rate at which the TLS releases quanta into 
the bath. We will use Eq. ([T5|) to generate simple ana- 
lytical results below. 

It can be seen that Eqs. pII)) -([T ^ do not constitute a 
closed set of equations. While the correlations in those 
equations are relevant to the full quantum solution of the 
problem, they turn out to play a negligible role in two 
cases : a) when the oscillator excitation is very low, since 
the TLS is then mostly in its ground state which then 
factors from the resonator state, and b) when the oscilla- 
tor excitation is very high, since in this case the TLS is in 
a fully mixed state. We therefore assume decorrelations 
such as {aa+) ~ {a){a+), etc., so that Eqs. pll)) - p^ now 
become 



(16) 



= 2ig{a){a++) - ig{a) - — (17) 
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Ti 
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which form a set of four closed equations if the conjugate 
of Eq. p7)) is included. The numerical solutions to these 
equations are substantially quicker to obtain compared 
to the master equation and will be discussed below with 
reference to their approximate analytic solutions. 



IV. EQUATIONS OF MOTION 
TREATMENT 



ANALYTICAL 



In order to identify the various regimes of loss and to 
organize the numerical solutions to the master equation 
[Eq. 0] , we now discuss the problem in terms of approxi- 
mate analytical solutions to the equations of motion. It is 
well known from previous studies that the JCM displays 
qualitatively different behavior in the regimes of weak 
[g < 1/Ti, I/T2) and strong {g > 1/Ti, I/T2) coupling^. 
We now consider the two regimes separately. 



A. Weak coupling : g < 1/Ti, I/T2 

For weak coupling the JCM dynamics are dominated 
by relaxation and dephasing processes, allowing for little 
coherence between the resonator and the TLS. Below we 
will see that in this case the TLS may be considered as 
being driven by the resonator at a time-dependent Rabi 
frequency 



n,^2g{n{t)y^^ 



(19) 



To find analytical solutions in this regime we consider the 
following argument. 

If the Rabi drive due to the resonator could be held 
constant, we could obtain steady state solutions for the 
TLS variables. However the drive depends on the number 
of resonator quanta, and thus itself changes with time. 
Nonetheless, if the photon number changes slowly com- 
pared to the rate at which the TLS relaxes to equilibrium, 
we can assume the drive to be quasistatic, evaluate the 
TLS variables in the steady state, and then finally solve 
for the resonator dynamics. The above scheme amounts 
to a Born-Oppenheimer approximation where the TLS 
dynamics arc considered as fast and the resonator dy- 
namics as slow. 

Solving Eqs. ([l6l)-(IT8l) for the TLS variables in the qua- 
sistatic state we find the population in the upper TLS 
state to be 



1 R^jt) 
2l + m{t)' 



where 



R{t) = ng^T\T2, 



(20) 



(21) 



is the Rabi frequency divided by the geometric mean of 
the longitudinal and transverse decay rates. Also 



■_m_ 



(22) 
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which shows that the coherence internal to the TLS is 
negligible at both small and large R{t). The quantity 
R{t), which is the ratio of the Rabi frequency to the geo- 
metric mean of the TLS decay rates, naturally delineates 
two qualitatively different types of behavior of the sys- 
tem, treated below. 



1. Unsaturated regime : R{t) <^ 1 

In this regime we find from the numerics that for ini- 
tially small Rabi frequency, i.e. i?(0) ^ 1, the resonator 
photon number decreases exponentially; thus it is true 
that R{t) < 1 for all time t. From Eq. (HOI), we find 



(23) 



Using this approximation wc can solve Eq. (jl5|) to obtain 
{n{t)) = (n(0))e-", (24) 
with the inverse of the effective decay rate given by 



(25) 



The first term on the RHS of Eq. (|25p corresponds to the 
time required to transfer a single quantum from the res- 
onator to the TLS. The second term corresponds to the 
time required for the TLS to release each quantum to 
the bath. The first time interval decreases with increas- 
ing coupling rate g and the dephasing time T2. However 
in the weak coupling regime this term is always larger 
than Ti. 

As shown in Fig.[T]the analytical expression of Eq. ([M)) 
agrees well with the numerical solutions of both the mas- 
ter equation and the Maxwell-Bloch equations. We note 
that exponential decay of oscillator energy corresponds 
classically to the motion of a pendulum damped by a 
viscous fluid such as air, and is often referred to as 'wet' 
friction. 

The analytical solution for the population (cr_|_+(t)) can 
be found self-consistcntly by using Eq. ([24| in Eq. ([23)) . 
The dynamics of the population are shown in Fig. [2] and 
exhibit the validity of the analytical solution except for 
the non-adiabatic behavior at early times when the TLS 
undergoes rapid excitation in our numerical simulations. 
It can be seen in this regime that the TLS is far from sat- 
uration and transports quanta efficiently into the reser- 
voir. 

We note that our R{t) <^ 1 theory agrees with the 
result for a qubit coupled to a TLS^. The expression for 
the qubit population relaxation rate for weak coupling 
in Ref. i.e. v\l^{, agrees with the first and dominant 
term of Eq. (^5)) with the appropriate substitutions = 
2g and j( = 1/Ti in the absence of pure dephasing {T2 — 
2Ti). 
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FIG. 1: This figure shows the time decay of the resonator plio- 
ton number from the initial value {n{0)) = 3 in the weak cou- 
pling, TLS-unsaturated regime for -R(O) ~ 0.2. The three sets 
of points correspond to the numerical solution of the master 
equation, the numerical solution of the decorrelated Maxwell- 
Bloch equations and the analytical solution of Eq. (|24p . re- 
spectively, and are in quite good mutual agreement. 
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FIG. 2: This figure shows the dynamics of the population in 
the TLS level |-f ) starting from the initial value {a++{0)) = 
in the weak coupling, TLS-unsaturated regime, for -R(O) ~ 
0.2. The numerical solution to the master and decorrelated 
Maxwell-Bloch equations match quite well, whereas the ana- 
lytical Born-Oppenheimer solution of Eq. (|23|l fails to capture 
the rapid initial excitation of the TLS. 



2. Saturated regime : R{t) S> 1 

In this regime we find from the numerics that for ini- 
tially large Rabi frequency, i.e. R{0) ^ 1, the resonator 
photon number decreases essentially linearly with time. 
For this regime we find from Eq. ()20[) . 



^ 1/2. (26) 
Using this approximation we can solve Eq. (|15|) to obtain 

1 



(nit)) = (n(0)) - —t. 



(27) 
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FIG. 3: This figure shows the time decay of photon number 
from the initial value {n{0)) = 500 in the weak coupling, TLS- 
saturated regime for -R(O) ~ 4.2. The three sets of points cor- 
respond to the numerical solution of the master equation, the 
numerical solution of the decorrelated Maxwell-Bloch equa- 
tions and the analytical solution of Eq. (|27|l . respectively. The 
first two agree well, while the analytical expression captures 
the behavior in the dominant linear regime. 
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FIG. 4: This figure shows the dynamics of the popula- 
tion in the TLS upper level starting from the initial value 
{(T++(0)) = in the weak coupling, TLS-saturated regime, 
for -R(O) ~ 4.2. The numerical solution to the master and 
decorrelated Maxwell-Bloch equations match quite well. The 
analytical Born-Oppenheimer solution of Eq. (|26p is constant 
and approximates the fact that the TLS is saturated most of 
the time. 



A plot of Eq. p7| shown in Fig. [3] matches well the nu- 
merical calculation in the regime of linear decay and gives 
correctly a slope of — l/2Ti. The linearity of the decay is 
violated only at very long times when the photon number 
becomes very low and R(t) is no longer large compared 
to unity. 

Wc note that linear decay of oscillator energy corre- 
sponds classically to the motion of an oscillator impeded 
by a frictional force always opposed to the oscillator's 
velocity but constant in magnitude, and is referred to 
as 'dry' friction. An analytically solvable example is of 
a mass on a spring, vibrating horizontally on a surface 
and damped by sliding friction proportional to the nor- 
mal forcG^ii^. Clearly the frictional force 'saturates' at 
a constant value in this instance. A comparison may be 
made between Fig. [3] and that for the energy decay for a 
classical oscillator under dry friction (Fig. 5 in Ref. HI), 
including the nonlinear decay towards the end. 

The population dynamics of the TLS are shown in 
Fig. m It can be seen that in this regime the TLS is 
almost always saturated and therefore can only accept 
quanta from the resonator limited by the rate at which 
it can release them into the bath. 



3. Loss tangent 

We now consider the loss tangent, i.e. the behavior of 
the resonator loss 1/Q where Q is the quality factor, as 
a function of (n(0)), the initial average photon number. 
Typically modelling dielectric loss requires including the 
loss model for many TLSs, due to a distribution in A 
arising from the amorphous nature of the material. This 
procedure is usually performed by averaging the solutions 



from the classical model and could be carried out numer- 
ically by extending the quantum model that we consider 
here. 

To begin we recall the classical result, which follows 
from Eq. 



1 



2U \l + n'^TiT2, 
with the classical Rabi frequency given by 
pcos OAqFo 



n 



2hE 



(28) 



(29) 



where Fq{^ Fq) is the amplitude of the classical electric 
field. In this model saturation of the loss occurs at the 
critical Rabi frequency given by fic = (TiT2)~^/^, which 
implies, via Eq. (P^ . a critical electric field 



Fr. 



2hE 



pcosOAoVTiT^' 



(30) 



For a proper comparison to the quantum results the 
classical field should be equated to the expectation value 
of the quantum field operator [see Eq. (|4])] in the coherent 
state I a) : 

^'o = («|Fo(« + a^)l") - 2(n(0))i/2j7^. (31) 
This yields, finally, 

1 2g^T2 



Qc w[l + i?2(0)]- 



(32) 



The classical loss tangent is usually plotted as a function 
of the dimensionless electric field 



^0 

Fc 



RiO). 



(33) 
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We note that the classical expression of Eq. (|32l) models a 
steady state measurement of the loss when the resonator 
is driven so as to ensure {n{t)) = {n{0)). 

Now we turn to the results of the quantum theory, 
keeping in mind that there is no drive in our theory. 
When the TLS is unsaturated, the photon number de- 
cay is exponential, and the loss tangent may be defined 
as (n) = —Lu(n)/Q, or 



{n{t)) 



{hit)) 



(34) 



where we have evaluated the quality factor at {n(t = 0)) 
in order to make a valid comparison to the classical case. 
This implies from Eq. ([24]) that 
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For very weak coupling, Q/j^^ — 2g'^T2/uj, which repro- 
duces the classical result of Eq. with i?(0) <C 1. 

When the TLS is saturated the photon number decay 
is not exponential. In this case we find that the overall 
decay is typically dominated by the linear regime, espe- 
cially for (n(0)) > 1 (see Fig.E]). Using Eq. 
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(36) 



Qfl,»i 2(n(0))Titj' 

which agrees with the classical result of Eq. ([32)) for 
i?(0) 3> 1, reproducing in particular the inverse scaling 
of loss with energy [oc (n(0))~"'^]. 

From the above analysis we may expect the onset of 
TLS saturation to occur at a critical number of quanta 
riw given by the intersection of the unsaturated [Eq. ([55])] 
and saturated [Eq. ([5S|) ] loss asymptotes. 
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2Tir 4,92^1 Ta ■ 



(37) 



where the second equality corrcponds to the very weak 
coupling regime and also to the condition R = 1 , in agree- 
ment with the classical theory. The quantity n^j denotes 
the number of resonator photons required to saturate the 
TLS; in the weak coupling regime we expect > 1. 

Loss tangents were numerically calculated using 
Eq. ([9]) for specific values of g, Ti and T2 and are shown in 
Fig. [5] (as a function of {n{0))) and Fig. [Hj (as a function 
ofFo/F,). 



B. Strong coupling : g > 1/Ti, I/T2 

In the strong coupling regime of the JCM, the coherent 
exchange of energy between the resonator and the TLS 
plays an important role, although the overall dynamics 
is still governed by damping and dephasing. 



FIG. 5: This figure shows loss tangents of the resonator as 
a function of initial photon number for various parameter 
regimes. The weak coupling quantum loss tangent (dashed 
curve) was generated using the values g = 0.2, Ti = 1,T2 = 
0.2. The strong coupling quantum and classical loss tangents 
(solid and dotted curves respectively) were both generated 
using the parameters g — 10, Ti = 1, T2 = 0.2. All curves are 
normalized to the weak coupling loss at low photon number 
given by Eq. (|35p . As expected, quantum theory implies a 
higher loss for strong coupling of the resonator to the TLS 
than for the weak. However the loss according to quantum 
theory is much less than predicted by the classical theory us- 
ing the strong coupling parameters. The photon number at 
which the loss saturates is n„(ns) for weak(strong) coupling 
and is shown as a dashed (solid) vertical line. 
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FIG. 6: This figure shows loss tangents of the resonator as 
a function of the dimensionless electric field of Eq. 1)33^ for 
various parameter regimes. The weak coupling quantum loss 
tangent (dashed curve) was generated using the values g = 
0.2, Ti = l,T2 = 0.2. The strong coupling quantum and 
classical loss tangents (solid and dotted curves respectively) 
were generated using the parameters g — 10, Ti = 1, T2 = 0.2. 
All curves are normalized to the weak coupling loss at low 
photon number given by Eq. (|35|) . A comparison to Fig. [S] 
shows that while saturation occurs at a lower photon number 
for strong coupling compared to weak, it occurs at a larger 
electric field. 
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1. Unsaturated regime : n <tins 

For small resonator photon numbers, we notice from 
the numerics that the coherent energy transfer between 
the resonator and the TLS occurs at about the vacuum 
Rabi frequency 2g > (1/Ti,l/T2). Thus the resonator 
photon number no longer changes slowly compared to 
the rate of TLS relaxation and the Born-Oppenheimer 
approximation employed above in the case of weak cou- 
pling ceases to be valid. 

However an analytic understanding may still be gained 
from Eq. ([9]) by observing from the numerics that for 
low photon numbers the system dynamics is limited to 
a small effective Hilbert space. The space is defined by 
only three states : |0, — ), |1, — ), and |0, +), where the first 
entry in each ket denotes the resonator number state and 
the second the TLS state^^. In this restricted manifold 
noting that (a^acr^) = — (a^a), we find from Eq. ^ 



— (aV_) = ig{{(J++) - (n)) - 



T2 



(38) 



Together with its conjugate, Eq. ([55)1 forms a closed set of 
equations with Eqs. (|T2|) and ((T4| . These equations can 
be solved subject to initial conditions set by the average 
photon number in the resonator (n(0)); the lack of TLS 
excitation {a++{0)) = 0, and the absence of correlations 
in the system {a''a^{0)) = 0. 

As a simple example let us first completely neglect pure 
dephasing (T2 = 2Ti). We find, from Eqs. ^ and 
the solution for the photon number. 



(39) 



= {n{0)) / cos^ gt, 
for the TLS population in the state 1+), 

{<j++{t)) = {n{0))e-'/^^^ sin' gt, 
and for the correlation 



(40) 



|(aV_(f))|2 = (n(0))2e-*/2Ti sin'2gt/4. (41) 

We note that these results agree with those for strong 
coupling between a qubit and a TLS^ as can be seen by 
using the substitutions v± = 2g,5uj = and ^{ = l/Ti 
in Eq. (4) of Reference d. 

In the presence of pure dephasing, Eqs. ([T2|) . p4| and 
([55t can still be solved analytically but the results are 
quite cumbersome and we do not present the formulas 
here. They are instead plotted in Fig. [7] and Fig. [5] and 
compare well with the numerical solutions of the full mas- 
ter and Maxwell-Bloch equations. When pure dephasing 
is included, our results do not agree exactly with the 
qubit-resonator predictions of reference H. The differ- 
ences can be traced to our use of the rotating wave ap- 
proximation (which discards terms from the equation of 
motion that oscillate at the bare TLS frequency lo but 
retains terms that oscillate with frequency 2g) and the 
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FIG. 7: This figure sliows tlie plioton number dynamics start- 
ing from the initial value (n(0)) = 0.005 in the strong cou- 
pling, TLS-unsaturated friction regime. The three sets of 
points correspond to the numerical solution of the master 
equation, the numerical solution of the decorrelated Maxwell- 
Bloch equations and the analytical solution to the system of 
equations given by Eqs. (|12|) . (|14|) and (|38p . respectively. 
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FIG. 8: This figure shows the dynamics of the population in 
the TLS level starting from the initial value {0-4-+ (0)) = in 
the strong coupling, TLS-unsaturated regime. The three sets 
of points correspond to the numerical solution of the master 
equation, the numerical solution of the decorrelated Maxwell- 
Bloch equations and the analytical solution to the system of 
equations given by Eqs. (fT2|) . (fl4ll and p8|) . respectively. 



use in Rcf. H of the secular approximation to the Bloch- 
Rcdfield equation (which also discards terms that oscil- 
late with frequency 2g). 

It may be noted that the results of this section were ob- 
tained from the master equation of Eq. ^ . That master 
equation was arrived at by first considering a TLS and 
its dissipation. A lossless oscillator was then coupled to 
the TLS. It has been pointed out that while this pro- 
cedure yields a good approximation in the case of weak 
coupling, a different master equation ought to be used 
for the case of strong coupling, one derived by first di- 
agonalizing the strongly coupled TLS and oscillator, and 
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then adding dissipation to the whole systeni^^"— . We 
have verified that the resonator loss for strong coupling 
calculated from such a master equation agrees with that 
presented in this article for the TLS-unsaturated regime. 



2. Saturated regime : n ^ rig 

We find from the numerics that for large photon num- 
bers, the resonator photon number decays linearly at a 
rate l/2T'i and the TLS is saturated. The resonator and 
TLS dynamics are similar to the case of weak coupling 
as shown in Figs. |3] and |4] respectively, and are well de- 
scribed by Eq. (|27|) and Eq. (|26|) respectively. 

However we note two aspects with respect to which the 
strong coupling dynamics differs from that in the weak 
coupling case. First, superposed on the overall decay of 
photon number and saturation of the TLS are oscillations 
of small amplitude which occur at early times as coherent 
transients, and later as revivals familiar to the JCM^. 
Second, saturation of the TLS does not occur at R{t) ~ 
1, in contrast to the case of weak coupling as discussed 
below. 



critical photon number given by the intersection of 
the unsaturated [Eq. (|42|) ] and saturated [Eq. (p6| ] loss 
asymptotes. 



3 / Ti 



(43) 



We note that Eq. (|43p is valid for very strong coupling, 
g ^ l/Tijl/Tlj. This explains the absence of g from 
Eq. (|43|) : deep in the strong coupling regime g is respon- 
sible for the reversible dynamics, while loss arises from 
irreversible processes. This argument also implies that 
the dimensionless quantity Us can then only depend on 
the ratio T1/T2 as verified by the final formula; for a large 
Ti the resonator easily saturates the TLS and a small T2 
implies that the coherent oscillations via which the TLS 
returns quanta to the resonator are quickly dcphascd. In 
the absence of pure dephasing {T2 = 2Ti) saturation be- 
gins to turn on at ris = 1; for T2 < 2Ti, n,, < 1. However, 
since it is difficult to transiently saturate the TLS with 
(ri(0)) ^ 1, Us docs not take arbitrarily low values. A 
numerical search yields the bound 



> 1/2. 



(44) 



3. Loss tangent 



CONCLUSION 



Although the qualitative shape of the loss tangent is 
similar to that of weak coupling, as can be seen from 
Fig. [5] and Fig. [6] there are crucial physical differences 
between the two. First, the asymptotic value for strong 
coupling when the TLS is unsaturated is found from the 
analytic solution of Eqs. ([12]), ([14]) and ([38|) to be 



1 












3uj ^ 






4( 



1 



1 

3f: 



Our prediction, Eq. indicates that for large g, the 

classical theory overestimates the loss in the unsaturated 
regime. Furthermore the result of Eq. (|42p can be used 
to show straightforwardly that the quantum theory pre- 
dicts a much greater loss in the unsaturated regime for 
large g (strong coupling) than for small g (weak cou- 
pling). Physically, this is reasonable - loss is higher for 
the resonator when it is strongly coupled to the TLS. 
Comparing Eq. to Eq. ([25)1 we see that the pure de- 
phasing plays opposite roles in the two cases : loss in the 
weak coupling case decreases with pure dephasing while 
it increases for strong coupling. In the saturated-regime 
all theories coincide, and follow Eq. which holds 

for strong as well as weak coupling. The loss curves are 
shown as a function of the dimensionless electric field of 
Eq. ([33|) in Fig. HI from which it can be seen that although 
the saturation photon number for strong coupling is lower 
than for weak coupling (Fig. [5]), the critical electric field 
at which saturation occurs is higher (Fig. |6]). 

Lastly, from the above considerations the crossover be- 
tween the two regimes can be estimated to occur at a 



Wc have quantized the classical model for dielectric 
loss in superconducting resonators due to two level charge 
fiuctuators, arriving thus at the Jaynes-Cummings model 
of quantum optics. For both the strong and weak cou- 
pling scenarios of this model we have identified regimes 
corresponding to nonsaturation and saturation respec- 
tively of the TLS. We have found that the quantum the- 
ory agrees with the classical in the regime of weak cou- 
pling. However for strong coupling, we find the quantum 
theory prescribes a substantially greater loss than the 
classical. Moreover, the photon number at which TLS 
saturation occurs is independent of the resonator- TLS 
coupling in the strong coupling quantum theory, in com- 
plete contrast to the classical theory. 

The numerical and analytical results presented in this 
paper can be used as a starting point for including further 
parameters in the loss model in order to make it more 
realistic, such as an external drive for the resonator, in- 
trinsic resonator loss, and reservoirs at non-zero temper- 
ature. In addition, we have in this work only considered 
the situation where the resonator is resonant with the 
TLS. The case of non-zero detuning is of interest since ex- 
periments typically include a distribution of many TLSs 
with varying detunings. Usually distributions in the vari- 
ables of orientation (9), energy (A) and tunneling bar- 
riers (Aq) also need to be considered in order to reach 
agreement with experiments^. Clearly there can be sev- 
eral additional regimes of behavior of the quantum loss 
model, depending on the relative values of the parameters 
introduced. While working in the time domain serves to 
illustrate the dynamics clearly, working in the frequency 
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domain can yield additional insights, and is relevant to future. We would like to thank F. Wellstood, M. Stouti- 
some experiments; such a project is envisioned for the more, and M. Khalil for discussions. 
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